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Abstract 

Several surfactant molecules self-assemble in solution to form long, cylin- 
drical, flexible wormlike micelles. These micelles can be entangled with each 
other leading to viscoelastic phases. The rheological properties of such phases 
are very interesting and have been the subject of a large number of experi- 
mental and theoretical studies in recent years. We will report our recent work 
on the macrorheology, microrheology and nonlinear flow behaviour of dilute 
aqueous solutions of a surfactant CTAT (Cetyltrimethylammonium Tosilate). 
This system forms elongated micelles and exhibits strong viscoelasticity at low 
concentrations (~ 0.9 wt%) without the addition of electrolytes. Microrheol- 
ogy measurements of G{ijj) have been done using diffusing wave spectroscopy 
which will be compared with the conventional frequency sweep measurements 
done using a cone and plate rheometer. The second part of the paper deals 
with the nonlinear rheology where the measured shear stress cr is a nonmono- 
tonic function of the shear rate 7. In stress-controlled experiments, the shear 
stress shows a plateau for 7 larger than some critical strain rate, similar to 
the earlier reports on CPyCl/NaSal system. Gates et al have proposed that 
the plateau is a signature of mechanical instability in the form of shear bands. 
We have carried out extensive experiments under controlled strain rate con- 
ditions, to study the time-dependence of shear stress. The measured time 
series of shear stress has been analysed in terms of correlation integrals and 
Lyapunov exponents to show unambiguously that the behaviour is typical of 
low-dimensional dynamical systems. 
PACS: 83.50.Fc, 83.50.Gd, 05.45.-a 

I. INTRODUCTION 



Rheology is the study of the deformation and flow of matter. Solids and fluids exhibit dif- 
ferent flow behaviours under shear. Solids store mechanical energy and are elastic, whereas 
fluids dissipate energy and are viscous. Complex fluids (e.g. colloids, polymers), owing 
to their larger length scales which result in low mechanical susceptibilities (~ 10^ Pa as 
compared to ~ 10^^ Pa in atomic systems), show very complex flow behaviour and are vis- 
coelastic. The relative proportion of elastic and viscous responses depends on the frequency 
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of the applied stress. For example, for entangled polymer solutions, the stress relaxation 
predominantly occurs by reptation dynamics with time scales r^ep- This system will be more 
elastic for uj > whereas it will be more viscous for uj < t~^. One should also note 
that for soft condensed matter, the elastic modulus under shear stress is much smaller than 
the elastic modulus under compressive stress whereas these two moduli are nearly equal for 
conventional atomic systems. 

Experiments on the rheology of matter involve the measurement and prediction of its 
flow behaviour. The method involves the application of a known strain or strain rate to a 
sample and the subsequent measurement of the stress induced in the sample or vice versa. 
The response of a viscoelastic material to an applied stress may be characterized as linear 
or nonlinear depending on the magnitude of the applied stress/strain rate. 



A. Linear Rheology 

The response of a material is linear when very small stresses (i.e. small compared to the 
spontaneous thermal fluctuations in the material) arc applied. If a small step strain 7 is 
applied to a deformable material at time t — ti, the stress induced in the material is given 
by o-.j.y(t) = G(t,ti)j{t) (Fig. 1), where G(t,ti) is the stress relaxation function. Here x 
and y define the velocity and velocity gradient directions, respectively. To linear order in 7, 
all other components of stress like a^x arc zero. Invoking time-translational symmetry, G 
depends on the time difference between t and ti i.e. G{t, ti) = G{t — ti). Exceptions to time- 
translational symmetry are found to occur in glassy systems which show aging behaviour. 
For an arbitrarily small applied strain rate 'y{t), the stress cr^ylt) induced in the material is 
defined as 

<J,y{t)^ f_^G{t-t'm')dt', (1) 

where G(t) is known as the memory kernel for shear response, 
(a) Oscillatory flow 

If a strain 7(t) = 706*"^* is apphed to the material, the stress developed in the sample 
can be out of phase with the strain by a phase angle S (Fig. 2). For a viscoelastic material, 
S lies between the limits S = (for a Hookean solid) and 5 — t:/2 (for a Newtonian fluid). 
Using Equation 1, the resulting stress may be written as 

(J.y{t)=lo\' G{t-t')iue'^''dt'. (2) 
/ —00 

Putting T = t-t', it can be written as a^y{t) = 7oe''^*G*(a;), G*{u;) = G'{u;) + iG"{uj), is 
the complex shear modulus given by 

G*{uj) = iuj\'^G{t)e-''^Ht. (3) 

The real part of G*{uj)^ called the storage modulus and denoted by G'{(jj), gives the elastic 
response of the material to the applied strain. It is the ratio of the stress in phase with the 
applied strain to the strain. The viscous response, deflned by the loss modulus G"{ijj), is 
the ratio of the out-of-phase stress component to the strain. Equation 3 implies that G'{ijj) 
is an even function whereas G"{ijj) is an odd function of uj. 
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(b) Steady shear 

When 7(t) = constant, the stress induced in the sample is given by a^y = 7 J G{t — 

t')dt' = rjo'j, where r/o is the zero shear viscosity given by r/o = J ^ G(t)dt. A non-zero value 
of rjo implies the presence of liquid-like or glassy dynamics. 

The simplest form of the response function G(t) is given by the Maxwell model, namely, 
G{t) = GoC"*/'^", or G*{lu) = GoIcotm/ {'^ + i^jtm), where Go is the elastic modulus and tm 
is the Maxwell relaxation time related to zero frequency shear viscosity rjo = GqTm- For a 

Newtonian fluid. Go — ^ oo, tm — ^ such that rjo remains constant. More generally, 

-t 

when there are many relaxation times present in the system, G{t) = J2j Gj e . For a 
continuous distribution of relaxation times, specified in terms of -P(t), 

G{t) = Goj^PiT)e~'/^dT = Goii{t). (4) 

For entangled polymer solutions, the stress relaxation occurs by reptation dynamics with 

time scale r^ep, for which = Y.n odd ^--p 

In a linear creep experiment, a step response is applied {a^y = for t < and a^y = o"o 
for t > 0) and strain 7(t) is measured, which is a solution of 

<yo= \\G{t-t'm')dt'. (5) 

The linear response of a viscoelastic material to an applied stress may be determined by 
a conventional frequency sweep experiment done using a rotating disc or concentric cylinders 
rheometer. The rheometer applies very small oscillatory stresses and calculates the resultant 
strain in the sample over the desired frequency range. For an applied angular frequency a;, 
the response that is in phase with the applied stress is used to calculate the storage modulus 
G'{uj) while the out of phase component gives the loss modulus G"{uj). 



B. Nonlinear Rheology 

Nonlinear rheology describes the response of a material to much larger stresses. As the 
name suggests, the strain induced in a sample varies nonlinearly with the applied stress in 
this regime. The nonlinear behaviour of a viscoelastic material in steady flow experiments 
is characterized by shear thinning or thickening, the presence of non-zero yield stress and 
normal stress differences, flow-induced phase transitions and the phenomenon of shear band- 
ing as shown schematically in Fig. 3. In nonlinear step strain experiments, if a large 
enough step strain 70 is applied to a sample, then the stress induced in the sample may be 
expressed by cr^yit) = joGniit — t',7o). The normal stresses under these conditions are no 
longer negligible and must be taken into consideration. For 70 — > 0, Gni G(t) as measured 
in linear rheology. 

Systems of giant wormlike micelles formed in certain surfactant solutions are known to 
show very unusual nonlinear rheology In steady shear, the shear stress saturates to a 
constant value above a critical strain rate 7 (as shown in Fig. 3 (d)) while the first normal 
stress difference increases roughly linearly with shear rate 0. Such behaviour is a signature 
of mechanical instability of the shear banding type and may be understood in terms 



3 



of the reptation-reaction model which involves the reversible breakage and recombination 
of wormlike micelles along with repation dynamics known for polymer solutions. Alterna- 
tively, the non-monotonicity of the flow curve has been attributed to the coexistence of two 
thermodynamically stable phases (isotropic and nematic) in the sheared solution 0. 

The flow curve may be measured under conditions of controlled stress or strain rate, 
and depending on the time interval between the collection of data points, we can obtain 
metastable or steady-state branches, respectively. In stress relaxation experiments, a con- 
stant step strain rate is applied to the sample in the nonlinear regime, following which the 
relaxation of stress in the sample is measured as a function of time. Alternatively, stress 
relaxation may be studied after cessation of a controlled strain rate that had been applied 
to the sample for a known duration. 



II. OUR EXPERIMENTS ON CTAT AQUEOUS SOLUTIONS 

We have studied the linear and nonlinear rheology of dilute aqueous solutions of the 
surfactant system CTAT (Cetyltrimethylammonium Tosilate) at 25°C. Above concentrations 
of 0.04 wt.% , and temperatures of 23°C 0, CTAT self-assembles to form cylindrical worm- 
like micelles which get entangled at concentration > 0.9 wt.%. The lengths of these wormlike 
micelles depend on the concentrations of the surfactant and the added salt, the temperature 
and the energy of scission of the micelle. The energy of scission is the excess free energy 
of a pair of hemispherical end caps relative to the rod like region containing an equivalent 
number of surfactant molecules. The number density of the elongated micelles of length L 
is given by 

Co{L) ~ ^exp(--^) (6) 

^ J-'avg 

where L is expressed in monomer units and 

L ~ 0O-5ea;p(4^) (7) 

where is the surfactant volume fraction and Esds is the energy of scission of the micelle. 
In these systems, stress relaxation occurs by reptation with time scale r^ep (the curvilinear 
diffusion of the micelle through an imaginary tube segment) as for conventional polymers 
and by the reversible scission (breakdown and recombination of micelles with time scale 
P]. The time scales r^ep and may or may not be comparable and depend on the surfactant 
concentration, presence of counterions in the solution and temperature, 
(a) Macrorheology measurements 

The frequency response of a viscoelastic material may be measured using a rheometer 
which consists of a device that can simultaneously apply a torque and measure the resultant 
strain. The in-phase and out-of-phase responses of the material to the torque are measured to 
calculate its elastic and viscous moduli, respectively. The instrument used by us is Rheolyst 
AR-IOOON (T.A. Instruments, U.K.) stress controlled rheometer with temperature control 
and software for strain rate control to measure the elastic and viscous responses of 1 wt.% 
CTAT between the angular frequency range of 0.03 rad/sec and 10 rad/sec. The rheometer 
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used was equipped with four strain gauge transducers capable of measuring the normal force 
with an accuracy of 10~^ N. The measurements were made using a cone-and-plate geometry 
of cone diameter 4 cm and angle 1°59". 

The linear regime of CTAT was first ascertained by looking for a range of stress values 
where the magnitude of the response functions were found to be independent of the applied 
oscillatory stress. The elastic and viscous moduli and the viscosity of CTAT 1 wt.% at 25°C 
were found to be constant for stresses between 0.05 and 0.1 Pa, oscillating at a frequency of 
0.1 Hz. Hence for the linear response measurement, 0.08 Pa was chosen as the amplitude of 
the oscillatory stress, oscillating between angular frequencies 0.03 rad/sec and 10 rad/sec. 
At frequencies higher than 10 rad/sec, the waveform depicting the strain becomes distorted, 
possibly due to the slip between the sample and the plates. This, therefore, limits the 
measurements till 10 rad/sec. Linear response measurements (Fig. 4 ) show that at the 
lowest frequencies, CTAT behaves like a viscous material, whereas in higher frequency runs, 
the behaviour is found to be predominantly elastic. The crossover is found to occur at 0.45 
rad/sec which corresponds to a relaxation time of 2.2 seconds. Cates et al have shown 
that for a system of wormlike micelles, like CPyCl/NaCl, G'{(jj) and G"{uj) are given by 
Maxwell model: 



G'M = G,^V//(1 + ^V), 
G"{uj) = G,ujtr/{1 + ujW). 



where tr = (ThTrep) ' ■ Figs 5 (a) and (b) show the least square fits of the data to the Maxwell 
model giving Go = 2.1 Pa and th = 2.2 sec. We find that for CTAT at concentration 1 wt.%, 
the fit is very poor. Further, the Cole-Cole plot (Fig 5 (c)) corresponding to the above data 
shows a deviation from the semi-circular behaviour expected in Maxwellian systems and 
shows an upturn at high frequencies. This deviation from Maxwellian behaviour is possibly 
due to the comparable values of r^ep and r?, in this system unlike in other wormlike micellar 
systems where the differences in the time scales {ti, << Tj-ep), lead to a 'motional averaging' 
effect. We have tried the Doi-Edward model where P] 

G'ico) = GoY: i^rn/pf/il + (coTn/pf) 

p odd 

G"{u) = G,Y: {^rr>/p)/{l + [ujTo/pf] 

p odd 

Here also, the fit with p = 1 and 3 is poor as shown in Fig. 6. The Doi-Edward model gives 
Go ~ 3 Pa and td ~ 1 sec. We also find that the Hess model which is given by 

G\uj) = ((r/o - r/oo)/re)^V,V(l + ^ V) 
G"{uj) = ((r^o - r7oo)/re)a;r,/(l + uj\') 

does not fit our data over the entire frequency range. In this model, rj^o is the high frequency 
shear viscosity and is a characteristic relaxation time of the system. It is likely that G'{uj) 
and G"{uj) can be fitted with a model calculated in low and high frequency regions separately. 



as is usually done in polymer literature 110]. Linear response measurements were also made 



for CTAT 1.9wt.% and CTAT 5wt.%. Interestingly, CTAT 1.9wt.% shows an anomalously 
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large relaxation time, whereas CTAT 5 wt.% is found to exhibit Maxwellian behaviour, as 
also seen in recent studies 0. 

(b) Microrheology measurements 

In recent years, micro-rheological techniques have been developed in addition to macro- 
scopic rheometry measurements using rotating disc or concentric cylinder rheometer. The 
basic idea behind microrheology is the tracking or manipulation of sub-micrometer particles 
immersed in the viscoelastic medium to be studied. Using magnetic beads as the probe 
particles, which can be manipulated by magnetic field gradients, the viscoelastic properties 
of F-actin networks |TT|] and the vitreous body of the eye \V2] have been measured. It is also 



possible to do microrheological measurements by a quantitative measurement of the mean 
square displacement < Ar^(t) > of the probe particles caused due to thermal fluctuations. 
This can be done either using laser interferometry with a resolution less than 1 nm |l3|jn 
or by diffusing wave spectroscopy flSl . 



The motion of the probe particle of radius 'a' may be described by a generalized Langevin 
equation given by : 

mv{t) = Ut)- j^at-T)v{T)dT (8) 

where mv{t) is the inertia of the particle, fnit) is the contribution due to electrostatic 
and Brownian forces on the particle. C, defines a time-dependent memory function which 
contributes to the viscous damping of the particle in the viscoelastic medium. The memory 
function C,{t) and are related by the following temporal autocorrelation function: 

< fR{0)fR{t) >= kBTCit) (9) 

where fc^ is the Boltzmann's constant and T is the temperature. 

In the frequency domain, the viscosity of the medium may be related to the frequency 
dependent memory function by the generalized Stokes-Einstein relation 

V{s) = 1^ (10) 

where s is the complex frequency given by s = ito. The complex viscoelastic modulus is 
given by |jl5 



G(s) = sri(s) = \— , , ■ ms] (11) 

The last term on the right hand side of Equation 11 is the contribution due to the inertia of 
the particle and can be omitted except at very high frequencies. < Ar'^{t) > of the probe 
particle is obtained from the intensity autocorrelation function which can be measured in 
diffusing wave spectroscopy experiments in the transmission or backscattering geometries. 
The Laplace transform of < Ar^(t) > is used to calculate G{s) using Equation 11. G{s) is 
then fitted to a functional form in s, which may then be used to calculate G*{u!), using the 
method of analytic continuation. 

We have used microrheology to estimate the G'^uj) and G"{uj) of an aqueous solution 
of CTAT of weight fraction 1%. The probe particles used are polystyrene colloidal parti- 
cles of diameter 0.23 /im dispersed in water at = 1%. Diffusing wave spectroscopy was 



6 



performed on the equilibrated sample using our light scattering setup consisting of a Kr^ 
ion laser (model 2020, Spectra Physics, U.S.A., excitation wavelength used 647.1 nm), a 
homemade spectrometer, photomultiplier tube (model R943-02, Hamamatsu, Japan), single 
photon amplifier discriminator (SPEX) and a MALVERN 7132 CE 64 channel correlator 
(Fig. 7). The light scattered Is(t) by the probe particles in the backscattering direction 
at a temperature of 25°C is used to measure the normalized intensity autocorrelation func- 
tion g2{t) =< Is{0)Isit) > / < /s(0) >^, as shown in Fig. 8 (a). For backscattering 
geometry, g2it) ~ Q-r(k'^<Ar^(t)>)^^^ which is used to get < Ar^(t) > as shown in 



Fig 8 (b), using F = 2. The parameter F is a constant depending on the polarisation of 
the scattered light and varies inversely with the transport mean free path V of the dif- 
fusing photon. Fig 8 (c) shows G(s) calculated using Equation 11, which was fitted to 
G{s) = Po+ pis-^-^^ + p2S^-'-^ + p3S°-^ +P4S |5[. Putting s = iuj, G'{uj) and G"{uj) of the 
dispersing gel are calculated as shown in Fig. 8(d). Comparison of Figs. 4 and 8 (d) shows 
similar magnitudes of the viscoelastic response functions obtained by macrorheology and 
microrheology methods. Further, in Fig. 8(d), the crossover of the viscous and elastic mod- 
uli occur at Uco ~0.4 rad/sec, indicating a relaxation time tr ~ 2.5 seconds for 1% CTAT 
at 25°C, similar to macrorheology measurements. It may be noted that microrheology may 
be used to calculate the frequency response of CTAT to much higher frequencies than the 
conventional rheometer experiment. 



III. NONLINEAR RHEOLOGY OF CTAT 

To study the nonlinear rheology of CTAT, we have measured the flow curve of CTAT 
1.35wt.% at 25° C as shown in Fig. 9 . The measurements are done under conditions of 
controlled stress. The data points are collected at intervals of 1 second, a value comparable 
to the relaxation time of the sample. The resultant branch of the measured flow curve is 
metastable, and its existence was demonstrated by Grand et al 0. The flow curve is found to 
saturate to a constant stress value above a critical shear rate 7c, while the flrst normal stress 
difference is found to increase linearly with shear rate. The plateau of the shear stress at 
high shear rates in CPyCl/NaSal has been interpreted by Grand et al as a characteristic 
feature of the flow curves of complex fluids that gives rise to a mechanical instability of the 
nature of shear banding . Shear banding results in the formation of bands of high and low 
viscosities in the sample, supporting low and high shear rates, respectively. However, the 
same phenomenon observed in CTAB/NaSal at a higher concentration has been explained 
by Berret et al as due to the coexistence of isotropic and nematic phases in the sheared 
sample. 

In addition to the measurement of the flow curve for our system, we have studied the 
stress relaxation in the sample after subjecting it to a step strain rate. At 25°C, on applying 
controlled shear rates whose values lie in the plateau region of the flow curve, the stress, 
instead of decaying to a steady state, is found to oscillate in time. Fig. 10 shows the time 
dependent stress relaxation in the 1.35wt.% CTAT sample at 25°C, on subjecting the sample 
to a step strain rate of 100s~^. The Fourier spectra of these time-dependent signals show 
time scales of the order of a few tens of seconds, which are an order of magnitude larger 
than tr. 

We identify the observed time-dependent behaviour as a manifestation of the mechanical 
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instability due to the formation of shear bands p|,|5|. Prehminary analysis of the time series 
obtained from the stress relaxation experiments done in the nonlinear regime shows the 
existence of positive Lyapunov exponents and finite correlation dimensions ll^(> 2 
at shear rates > 75 s~^), which points to the existence of deterministic chaos in sheared 
aqueous solutions of CTAT. The Lyapunov exponent characterizes the divergence of stress 
trajectories in the system, whereas the correlation dimension gives us information about the 
geometry of the attractor on which the trajectories in phase space asymptotically lie. On 
increasing the temperature of the sample to 35°C, and on maintaining the same shear rates 
as in the previous experiments done at 25°C, the time dependent oscillations in the stress 
relaxation are found to disappear completely. This is in accordance with previous studies on 
the temperature dependence of the flow curve of CPyCl/NaSal which shows a decrease 
in the width of the plateau with increasing temperature. The disappearance of the time- 
dependent behaviour in sheared CTAT at higher temperatures is thus a direct consequence 
of the disappearance of the shear bands in the sample . We have done extensive studies on 
the time-dependence of the stress relaxation of dilute, aqueous, sheared solutions of CTAT 
by doing more elaborate analysis of the time-series obtained from our experiments . 
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FIGURES 



FIG. 1. (a) shows step strain applied at time ti, with a magnitude small enough to lie in the 
linear regime, (b) shows the corresponding relaxation of stress. 

FIG. 2. 7 is the shear applied to a viscoelastic material, a is the resultant stress in the 
material, delayed by a phase angle 6. The complex shear modulus may be then written as 
G*{lo) = G'{u;)+iG"{io), where G'{uj) is the in-phase elastic modulus and G"{uj) is the out-of-phase 
viscous modulus. 

FIG. 3. (a) shows Newtonian {axy = rijo), shear thinning and shear thickening behaviours seen 
in viscoelastic fluids, (b) shows the flow curve of a Herschel-Bulkley plastic with a flnite yield stress 
(Jq. At axy < (To, the material behaves like a solid. At > ao, axy = a^ + Kp'j'^, which rcdiices 
to the Bingham equation for p=l. (c) shows the flow curve of a system undergoing a flow- induced 
phase transition, characterised by a sudden jump in a as shown by the dotted line, (d) shows a 
flow curve showing a plateau region, which is a signature of shear banding. 

FIG. 4. Measurement of the elastic modulus G'{u!) and the viscous modulus G"{uj) of lwt.% 
CTAT by the method of macrorheology, using a cone and plate rheometer, with applied stress = 
0.08 Pa at 25°C. 

FIG. 5. (a) G'{io) and (b) G"{io) of lwt.% CTAT using a rheometer (data same as shown in 
Fig. 4). The solid lines show the least square fits to the Maxwell model, (c) shows the Cole-Cole 
plot which deviates from the semicircle exhibited by Maxwellian systems. 

FIG. 6. Data same as in Fig. 4 and the corresponding fits to the Doi-Edward model shown 
by solid lines. 

FIG. 7. Our light scattering setup in the transmission geometry : PI and P2 are the analyser 
and polariscr respectively, LI and L2 arc convex lenses of f = 20 and 30 cm respectively, P.H. is 
a pin hole, S the sample that scatters light and Inv is an inverter circuit. Ml and M2 are plane 
mirrors that steer the incident laser beam. 

FIG. 8. The microrheology results: (a) shows the correlation function in the backscattering 
direction, (b) the < Ar^(t) >, (c) the complex modulus G(s) and (d) shows the calculated values 
ofG'{u) and G"{u). 

FIG. 9. The metastable branch of the flow curve of 1.35wt.% CTAT, measured under condi- 
tions of controlled stress . 

FIG. 10. The time-dependent relaxation of stress in 1.35wt.% CTAT, on subjecting the sample 
to a constant step-strain rate of 100 s~^. 



10 









(a) 


Y 










time 

(b) 


c 









time 

Fig. 1 

H 

Fig. 2 

11 





(a) / 

shear / y 
thickening/ y ^ 


(b) 




■ 




a 

xy 


>^ shear 






thinning 

ri 






(c) 


(d) 


a 

xy 


/ 

/ 





shear rate y 



Fig. 3 



10" - 



10' 



1 1 M M| 1 

r G(co) , 


1 1 1 M M 1 1 1 1 M M 


; + 




t 




1 1 


1 



10"' 10" 
CO (rad/sec) 



10 



Fig. 4 



12 




0) (rad/sec) w (rad/sec) 




Fig. 5 




(rad/sec) 



Fig. 6 



13 





Fig. 8 



14 




Fig. 9 




500 1000 1500 2000 2500 
time (sec) 



Fig. 10 



15 



